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1 INTRODUCTION 



One of the major problems in the analysis of diffraction data of 
disordered systems is the lack of any general method for develop- 
ing structural models that agree quantitatively with the experi- 
mental data. Most analyses are mainly qualitative and are based 
on a few data features, such as peak positions and coordination 
numbers derived from radial distribution functions. 

Sometimes the results of simulations are in good agreement with 
experiment (though comparison is made with the radial distri- 
bution functions rather than the structure factors), but usually 
the agreement is only qualitative, and occasionally there are ma- 
jor differences, especially in the descriptions of the liquid state 
of the 3(i transition metals. Usually, the disagreement between 
model and experimental results is due to the choice of interactive 
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potentials [1-6]. In this case, the inverse methods [7-12], in which 
the effective pair potentials are derived from the experimental 
structural data, is one of the best methods to compute models, 
structural properties of which quantitatively agree with the ex- 
perimental data. It is also well recognized that the effective pair 
potentials estimated from the experimental structural data are 
quite useful, because they are considered to include, more or less, 
the particular features of liquid of interest. 

The other problem of the simulation the disordered system is the 
analysis of the obtained information, which represent the com- 
plete list of coordinates and velocities of the atoms. Such infor- 
mation allows us to describe the realization of the current state 
of the system completely. One of the methods of solution of such 
problems is the statistical analysis of the Voronoy polyhedra built 
for each atom of the system [13,14]The analysis makes possible 
the description of the local order of the simulated system. 

In the present work the effective pair potentials for liquid 3d tran- 
sition metals Fe, Co and Ni were derived at the range of tem- 
peratures by means of the inverse method proposed in jHI from 
the experimental structural factors jIE|. To test the potentials, 
the dynamic properties, calculated by means of the molecular 
dynamics (MD) method, have been compared with their exper- 
imental values. Finally, the statistical analysis of Voronoy poly- 
hedra, built on the MD configurations of the investigated metals, 
has been carried out. 

2 METHOD 

2. 1 THE SELF- CONSISTENT METHOD 



The iterative self-consistent method originally proposed by Reatto 
et al. P is based on the computer simulations and a formally ex- 
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act result of the integral equations theory, which relates the effec- 
tive pair potential with the pair and direct correlation functions 
g{r) and c(r), respectively: 



where P = l/Zc^T, B{r) is so-called bridge-function. The func- 
tions g{r) and c(r) can be easily obtained from the experimental 
structural factor S^xpiq)- However, the 6r?(i(/e-function B{r) can 
not be obtained the same way. This is why in the original Reatto 
method the bridge-function of the hard sphere system Bhs{i^^ v)^ 
where the external parameter rj is the packing factor, was used as 
an initial estimation for the bridge-function, and potential (po{r) 
correspondently. As was pointed in the work [9j the convergence 
of the iterative method depends on the initial approximation of 
the potential <^o{r). To avoid the dependence on any external pa- 
rameters in the presented work the mean- field potential fMpif): 
which characterizes certain self-consistent field of particles, where 
the force acting on the particle of the system is the function of 
all particles arrangement, was used as the initial potential: 



The procedure of deriving effective potential, (p{r), is iterative. 
The equilibrium atomic configuration simulated by means of the 
molecular dynamics method with given potential (pi{r) corre- 
sponds to the pair correlative function gi{r), direct correlative 
functions, Ci{r) and structural factor, Si{q). The next iteration 
constructs the corrected potential: 

= 9exp{r) - 1 - Cexp{r) - In g^^p{r) + 5,(r), 
where the bridge-function Bi{r) exactly corresponds to the po- 
tential (fi{r): 



f](l)[r) = g{r) — c(r) — 1 — lng{r) + B{r) 



(1) 



^MF{r) 



ln{g{r)). 



Bi{r) = /3(pt{r) - gi{r) + ^(r) + \ng^{r). 
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The iterative process is repeated until the discrepancy x be- 
tween the criterion (experimental) and simulated functions be- 
comes smaller than the desired accuracy: 

1 1 M 

X=j^ Ufn'ir) - 9:{r)? + - EK^^(^) - S;^{q)]', (2) 

where and M are the numbers of points in the diagrams of the 
functions g{r) and S{q). 

2.2 COMPUTATIONAL DETAILS 



The effective potentials, obtained by means of the inverse method, 
reproduce the structural properties of the modelled system with 
good agreement with the experimental data. To check the ade- 
quacy of the simulated systems to a real ones it is also neces- 
sary to compare the dynamic properties of the system with their 
experimental values. Therefore, the molecular dynamics (MD) 
method, which enables the calculation of both structural and dy- 
namic properties of the system, was used in the work presented. 

The MD simulations of the 3(i liquid transition metals have been 
performed under microcanonical [NVT) ensemble. A^=5000 spher- 
ically symmetric particles, which interact by means of the effec- 
tive potential (^(r) , were arranged in a cubic box with the imposed 
periodical boundary conditions. The size of the box was deter- 
mined from the condition that the simulated system density is 
equal to the experimentally measured density of the alloys jTSj. 

The dynamic properties such as self-diffusion coefficient D and 
shear viscosity rj are known to be sensitive to the specified poten- 
tial of interaction in a system. To investigate the validity of the 
derived potentials the coefficients D and 77 have been calculated 
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by standard Green-Kubo (GK) formulas [r7|: 

^ = :tt7 / E < v^itoMt) > dt, (3) 



v = ^j < r^hw^t) > dt, (4) 

where a(3=xy, yz, zx and (jf^ is defined by follow expression: 

I j>i ar^j 

where Vi{t) is the velocity of i-ih particle with mass at the 
time t. 

2.3 STATISTICAL ANALYSIS OF THE VORONOY POLYHEDRA 



The sensitivity of the statistical analysis of Voronoy polyhedra to 
the local structure is particularly helpful for understanding the 
temperature evolution of the structure of disordered systems. The 
Voronoy polyhedron (VP) is defined as the sub-volume whose in- 
terior is closer to a specific atomic vertex than to any other vertex 
and may provide useful information about the local atomic en- 
vironment. The most widely applicable forms of the statistical 
analysis of the VP of disorder systems is the analysis of its met- 
ric and topological properties - distribution of VP volumes, face 
areas, topological indexes etc. However, it is more convenient to 
perform the analysis by means of the non-dimensional measure - 
sphericity coefficient Kgph'- 

i^sph = ^3 • (5) 

Here, V is the volume, and S is the surface area of the given 
polyhedron. Thus, the measure K^^^h is constructed to be unity 
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for regular sphere, and defines a deviation of the shape of the 
given VP from regular sphere. Since each VP is the geometrical 
image of the local environment of the given atom, the measure 
Kgph will characterize homogeneity of the its local environment. 
The statistics of the VPs, as the distribution of the sphericity co- 
efficients N{Ksph): will characterize whole system and for homo- 
geneous systems such distribution must have symmetrical Gauss- 
like shape. 



3 RESULTS AND DISCUSSIONS 



In the work the experimental structural factors Sexp{q) of fiquid 
iron at T = 1833, 1873, 1923, 2023, cobalt at 1823, 1873, 1923, 
2023 and nickel at 1773, 1873, 1923 and 2023 K [IB] were used 
as the criterion functions for inverse self-consistent method. The 
obtained effective pair potentials cp{r) have reproduced the struc- 
tural properties of the simulated systems with good agreement to 
the experimental ones. The deviation x between simulated and 
experimental functions reaches only 10^^ (fig. 1). 

The derived potentials, presented on fig. 2 with the experimental 
and simulated pair correlation functions, show long-ranged oscil- 
lations. For metallic systems such oscillatory behavior is inter- 
preted as the dielectric screening of ions by conducting electrons 
and can be associated with the Friedel's oscillations. 

There is a hump in the main minimum of the effective potentials, 
which becomes less pronounced with the temperature rise. At 
motion along the 3d row from Fe to Ni the intensity of the 
hump decreases as well. 

It is necessary to note, that position of minimum of the derived 
potentials rmin is shifted with respect to the position of the first 
maxima of RDFs vq. We believe, it can evidence a noticeable con- 
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tribution of the distant neighbors in the effective interaction. The 
ratio ro/rmin tends to be unity at high temperatures and at the 
motion along 3d row, that testify that the distant contributions 
weaken with temperature rise. 

To investigate the vahdity of the derived potentials in the present 
work the self-diffusion D and shear viscosity 77 coefficients were 
calculated for investigated melts. The obtained values are pre- 
sented on fig. 3 as the function of temperature. There is reason- 
able agreement between calculated and available experimental 
and theoretically predicted data [18-29]. 

Thus, the MD configurations, obtained using the effective poten- 
tials, can pretend on the realistic description of the structure of 
the liquid 3d transition metals, since reproduce both structural, 
and dynamic properties of the melts well. 

The distributions of the sphericity coefficient N{Ksph) were cal- 
culated Fe, Co and Ni at different temperatures and presented 
in fig. 4a, b, c. It is obvious that the shape of the curves N{Ksph) 
evolves with the temperature rise. Their maxima are shifted to- 
wards the low values of Kgph that testify the distortion of the local 
atomic environment in simulated liquids. At relatively low tem- 
peratures the shoulders on the curves were obtained. This fact 
can be evidence of the existence of some competitive types of 
the local order. However, with the temperature rise these shoul- 
ders disappeared and shape of the distributions N{Ksph) becomes 
more symmetrical, which says about more homogeneous structure 
of the liquid at this temperatures. 

To investigate the temperature infiuence on the local structure 
of the alloys it is convenient to determine the most probable 
values of the sphericity coefficient K^pj^, which characterizes the 
prevailing types of the local order in the structure and estabhsh 
the dependence of such values on the reduced temperature r: 



where and T are the melting and current temperatures, re- 
spectively. 

Fig. 4d demonstrates that at small values of r the maxima of 
the curves N{Ksph) for all simulated metals lie near the value 
Ksph = 0.7. We suppose, that this value corresponds to the com- 
mon symmetry of the local order, which is intrinsic to the dense 
liquids. 

The slope of the dependencies Kl^^ij) changes at r*, where the 
distributions N(Ksph) get more asymmetrical shape. It should 
note, that at r > r* the points of the dependencies for all in- 
vestigated metals can be approximated by the common linear 
law (dashed line). The temperature evolution of the distributions 
N(Ksph) ai T > T* allows us to suppose that the differences be- 
tween the competitive types of the local order in alloys vanish 
and the structure of liquid becomes more homogenous. 



4 CONCLUSIONS 



The results of the present work allow us to suppose, that the ob- 
tained effective potentials of liquid 3d transition metals at given 
temperatures can reproduce their real structure, since the model 
adequately reflect both structural, and dynamic properties of the 
melts. 

The performed statistic analysis of the Voronoy polyhedrons al- 
low us to estimate the temperature range of inhomogeneity — 
T* (where T* = Tmij* + 1) see fig. 4d). Within of this range 
metallic liquid contains several competitive types of the local or- 
der. Out of this range T > T* the structure of liquid becomes 
more homogenous and the character of the local order evolution 
becomes common for all the investigated metals. According to our 
estimation, the boundary of the temperature range of inhomo- 
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geneity is for iron T* = 1873 K{r* = 0.024), for cobalt T* = 1923 
K{t* = 0.061) and for nickel T* = 2023 K{t* = 0.114). We be- 
lieve, that the estimations can be useful on the choice of techno- 
logical parameters of supercooling of 3d metal-based alloys. 
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5 Figure captions 



Figure 1. The experimental (open circles) and simulated (solid 
lines) structural factors S{q) for liquid (a) Fe, (b) Co, (c) Ni 
at different temperatures. Deviation x between the functions 
reaches only 10~^. 

Figure 2. The experimental (open circles) and simulated (solid 
lines) pair correlation functions g{r) with the effective potentials 
(f{r) derived for (a) Fe {h)Co and (c) Ni from the experimental 
data pEj- The nearest neighbors distance ro is shifted with respect 
to the position of main minimum of the potential rmin- Since 
the ratio ro/rmin tends to be unity at the temperature rising 
and at the motion along 3d row from Fe to Ni, it means the 
contributions of the distant neighbors to the effective interaction 
weaken. 

Figure 3. Values of the self-diffusion coefficient D and shear vis- 
cosity r] (open squares) calculated for the MD configurations of 
Fe, Co and Ni using the Green-Kubo (GK) relations, and avail- 
able experimental and calculated data [1, 18-29]. 

Figure 4. Distributions of the sphericity coefficient N{Ksph) cal- 
culated for liquid (a) iron, (b) cobalt, (c) nickel in the investigated 
temperature range; (d) The most probable values of the spheric- 
ity coefficient i^f^/j as the function of the reduced temperature 
r. The slope of the dependencies Klp^ir) is changed at r*. The 
behavior of Kgpf^{r) at r > r* can be approximated by linear low 
(dashed line), which is common for all metals investigated in the 
present work. 
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